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ABSTRACT 

When embedded in dense cluster cores, intermediate mass black holes (IMBHs) acquire close stellar or 
stellar-remnant companions. These companions are not only gravitationally bound, they tend to hierarchically 
isolate from other cluster stars through series of multibody encounters. In this paper we study the demographics 
of IMBH companions in compact star clusters through direct A-body simulation. We study clusters initially 
composed of 10^ or 2 x 10^ stars with IMBHs of 75 and 150 solar masses, and follow their evolution for 6-10 
Gyr. A tight innermost binary pair of IMBH and stellar object rapidly forms. The IMBH has a companion 
with orbital semi-major axis at least three times tighter than the second-most bound object over 90% of the 
time. These companionships have typical periods of order years and are subject to cycles of exchange and 
destruction. The most frequently observed, long-lived pairings persist for ^10^ yr. The demographics of 
IMBH companions in clusters are diverse; they include both main sequence, giant stars, and stellar remnants. 
Companion objects may reveal the presence of an IMBH in a cluster in one of several ways. Most-bound 
companion stars routinely suffer grazing tidal interactions with the IMBH, offering a dynamical mechanism to 
produce repeated flaring episodes like those seen in the IMBH candidate HLX-1. Stellar winds of companion 
stars provide a minimum quiescent accretion rate for IMBHs, with implications for radio searches for IMBH 
accretion in globular clusters. Einally, gravitational wave inspirals of compact objects occur with promising 
frequency. 

Subject headings: black hole physics - stars: kinematics and dynamics - galaxies: star clusters: general - 
methods: numerical 


1. INTRODUCTION 

Globular clusters (GCs), dense young star clusters, and the 
compact nuclear star clusters of low-mass galaxies are envi¬ 
ronments of extreme stellar density. Their typical sizes of par¬ 
secs and masses > 10 ^Mq imply stellar number densities be¬ 
tween 10^ and 10^ times that of the solar neighborhood and 
velocity dispersions of tens of km s“^ (Heggie & Hut 2003; 
Neumayer & Walcher 2012). With similar ratios of stellar 
mass to system mass, stellar escape velocity to velocity dis¬ 
persion, and stellar lifetime to system relaxation time, these 
dense stellar systems have broad dynamical similarities de¬ 
spite their disparate environments. These dense clusters are 
environments of intense stellar interactivity, where single and 
binary stellar evolution and gravitational dynamics intertwine 
to shape the long-term evolution of the cluster as a whole 
(Heggie & Hut 2003; Aarseth 2003; Kalirai & Richer 2010). 

There has long been speculation that these dense stellar sys¬ 
tems might harbor intermediate mass black holes (IMBHs) 
with masses between that of the known populations of stellar 
mass and supermassive black holes (BHs) (e.g. Lightman & 
Shapiro 1978; Miller & Hamilton 2002). If the Mbr-ct rela¬ 
tion of supermassive BH mass and velocity dispersion were to 
extrapolate to objects as small as GCs and nuclear clusters, it 
would imply that these stellar systems should host BHs in the 
IMBH mass range (e.g. Eerrarese & Merritt 2000; Lutzgen- 
dorf et al. 2013c; Baldassare et al. 2015). In young clusters 
of low metallicity, evolving massive stars might preferentially 
leave behind massive remnant BHs, with masses > IO^Mq 
(M apelli et al. 2013c; Spera et al. 2015). Through dynamical 
interaction these BHs could acquire and accrete mass from 
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companion objects (Mapelli et al. 2013c; Mapelli & Zampieri 
2014; Ziosi et al. 2014) implying a self-consistent channel by 
which the BH could grow into the IMBH mass range as the 
cluster evolves. 

There is growing evidence that some nuclear star clusters 
in low-mass galaxies host BHs in the IMBH mass range (e.g. 
Neumayer & Walcher 2012). These sources have been dis¬ 
covered primarily in the fraction of systems where the central 
BH is active. Eor example, Seth et al. (2008) study the co¬ 
incidence of BH activity and nuclear star clusters, and And 
that a large fraction of all nuclear star clusters show signs 
of BH activity (> 50% of those with masses > IO^Mq). 
Seth et al. (2010) found a BH mass of ^ 5 x 10 ^Mq in the 
NGC 404 nucleus. By studying the active BH fraction in 
a large sample of Sloan Digital Sky Survey dwarf galaxies, 
Reines et al. (2013) And a ^ 0.5% active fraction in their 
sample of 10^-^ - IO^-^Mq dwarfs with median BH masses of 
^ 2 X IO^Mq. Reines et al. (2013) note that sensitivity limits 
their ability to detect active BHs of less than IO^Mq, even if 
they are emitting near their Eddington luminosities. Even so, 
their data suggest a similar active fraction among these low 
mass BHs to that of IO^Mq BHs. More recently, Baldassare 
et al. (2015) report a 5 x IO^Mq IMBH in the dwarf galaxy 
RGG 118 accreting at 1% its Eddington luminosity. The fact 
that this BH lies along the extrapolation of the Mbh “ cf rela¬ 
tion to low masses suggests that BHs in this mass range may 
be common in dwarf galaxies, with their detection hampered 
only by their low characteristic luminosities. 

The evidence for the presence or absence of IMBHs in GCs 
remains contested. An IMBH, unless illuminated through ac¬ 
cretion, would be be dark (e.g. Strader et al. 2012; Ramirez- 
Ruiz & Rosswog 2009; Volonteri et al. 2011). This leaves 
only the IMBH’s gravitational interaction with surrounding 
stars as a tracer of its presence. Searches for peaks in the 
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central velocity dispersion or surface brightness of GCs due 
to orbital motion around an unseen IMBH have returned con¬ 
tested results (e.g. Gebhardt et al. 2005; Lanzoni et al. 2007; 
Noyola et al. 2008; Miocchi 2010; Noyola et al. 2010; Noyola 
& Baumgardt 2011 ; Lutzgendorf et al. 2011; Jalali et al. 2012; 
Lutzgendorf et al. 2012a,b, 2013b; Feldmeier et al. 2013; Ka- 
mann et al. 2014; Lutzgendorf et al. 2015), with a fundamen¬ 
tal challenge being that unless the IMBH is a large fraction 
(e.g. > 1%) of the cluster mass very few stars will trace or¬ 
bits dominated by the IMBH’s gravity (e.g. Vesperini & Trend 
2010; Umbreit & Rasio 2013). This challenge has lead to the¬ 
oretical and observational efforts to identify other dynamical 
tracers of the presence of IMBHs in GC systems (e.g. Baum¬ 
gardt et al. 2005; Heggie et al. 2007; Hurley 2007; Gill et al. 
2008; Pasquato et al. 2009; Pasquato 2010; van der Marel & 
Anderson 2010; Beccari et al. 2010; Samra et al. 2012; Leigh 
et al. 2014). 

In this paper, we consider a stellar dynamical implication of 
the possible presence of IMBHs in dense cluster cores: stars 
bound tightly to the IMBH. IMBHs in clusters acquire a num¬ 
ber of gravitationally bound stars that make up their sphere 
of influence (Peebles 1972; Bahcall & Wolf 1976; Preto et al. 
2004; Baumgardt et al. 2004a,b). These stars orbit under the 
combined influence the IMBH and all of the other bound stars 
(Frank & Rees 1976). The stellar dynamics of the relatively 
small number of stars in the sphere of influence region is com¬ 
plex, and relies on the interface of stars’ orbital, dynamical re¬ 
laxation, and evolution timescales. In simulation, the most ac¬ 
curate (but also computationally expensive) approach to treat¬ 
ing the multiphysics sphere of influence dynamics is direct 
A-body evolution of the stars’ orbital motion (Aarseth 2003). 
Companion objects to IMBHs have been mentioned in previ¬ 
ous studies (e.g. Baumgardt et al. 2004b; Blecha et al. 2006; 
Baumgardt et al. 2006; Konstantinidis et al. 2013; Leigh et al. 
2014; Giersz et al. 2015), however, in this paper we present 
the first systematic study of IMBH companions in a dense 
cluster context that takes advantage of the accuracy of direct 
A-body integration. 

Stellar objects or stellar remnants tightly bound to an IMBH 
are of extreme astrophysical interest because of their ability to 
reveal the presence of the otherwise dark BH. The category of 
ultraluminous X-ray (ULX) sources is particularly interesting 
because their X-ray luminosities of > 10^^ erg s“^ preclude 
an Eddington limited stellar mass BH source (e.g. Humphrey 
et al. 2003; Pasham et al. 2014). In dense clusters, stars can be 
scattered into close, disruptive encounters with IMBHs, pow¬ 
ering luminous but brief flares of ULX luminosity (Ramirez- 
Ruiz & Rosswog 2009). If, instead, stars are captured into 
tight orbits and induced to Roche lobe overflow, persistent 
ULXs could conceivably be the very high mass analogs X-ray 
binaries (Hopman et al. 2004; Li 2004; Portegies Zwart et al. 
2004; Blecha et al. 2006; Patruno et al. 2006; Baumgardt et al. 
2006). If less tightly bound to the IMBH, a stellar companion 
could still fuel accretion through its stellar wind (Cuadra et al. 
2005, 2006, 2008; Miller et al. 2014). 

The source HLX-1 (Farrell et al. 2009) is particularly dra¬ 
matic, because its lightcurve shows repeated flaring episodes 
during which the source undergoes state transitions, strongly 
suggesting an IMBH as the central engine (Godet et al. 2009; 
Farrell et al. 2011; Davis et al. 2011; Servillat et al. 2011; 
Webb et al. 2012, 2014; Farrell et al. 2014). The roughly 
periodic flaring behavior has also lead to an association of 
the behavior with a IMBH plus star system in an eccentric 


orbit - where close pericenter passages might account for 
the brightening (Lasota et al. 2011; Soria 2013; Miller et al. 
2014; Godet et al. 2014). The association of the object with 
a compact stellar system strengthens the interpretation that 
the accretion flares might be powered by a stellar compan¬ 
ion (Wiersema et al. 2010; Farrell et al. 2012; Mapelli et al. 
2012, 2013b,a; Soria et al. 2013). An avenue for these re¬ 
peated, grazing tidal interactions with supermassive BH hosts 
has been suggested by MacLeod et al. (2013), where stars 
evolve to ‘spoon-feed’ portions of their envelope to the su¬ 
permassive BH over many passages. In the galactic center 
context, Guillochon et al. (2014) have proposed this channel 
as a dynamical mechanism to produce the tightly-bound G2 
gas cloud. However, the dynamical processes that lead an 
IMBH, like the putative BH in HLX-1, to acquire then strip 
mass from a tightly bound companion remain poorly under¬ 
stood. 

This paper studies the dynamical processes by which 
IMBHs in dense clusters acquire, retain, and lose close com¬ 
panion stars. We And that a stellar object among the sphere of 
influence members tends to segregate to substantially tighter 
separations than the other bound stars. These companion ob¬ 
jects form a hierarchically separated binary with the IMBH, 
and persist until they are replaced in an exchange interaction 
or destroyed through direct interaction with the IMBH. We 
characterize the orbital demographics and statistics of these 
companion objects using the A-body simulations presented in 
this paper. In Section 2 we describe our direct A-body sim¬ 
ulation method and parameter choices. In Section 3 we ex¬ 
plore the stellar dynamics, populations, and processes govern¬ 
ing tightly-segregated companions to the IMBH. In Section 4 
we examine the dependence of our simulated populations of 
companion on some key cluster parameters. In Section 5 we 
discuss some astrophysical implications of our results, par¬ 
ticularly the role that these tightly bound companions might 
play in revealing the presence of otherwise hidden IMBHs. 
Finally, in Section 6 we conclude. 

2. DIRECT iV-BODY SIMULATIONS OF CLUSTERS HOSTING 
IMBHS 

2.1. Simulation Method 

In this work we use the direct-summation A-body star clus¬ 
ter dynamics code NBODY6 (Aarseth 1999, 2003) in a ver¬ 
sion enabled for GPU-accelerated calculations (Nitadori & 
Aarseth 2012)^ to follow the evolution of star clusters in a 
realistic tidal held following the implementation of Heggie 
et al. (2006) and subsequent papers (Trent! et al. 2007b,a; 
Trent! 2008; Trent! et al. 2008; Gill et al. 2008; Pasquato et al. 
2009; Trent! et al. 2010; Vesperini & Trent! 2010; Trent! & 
van der Marel 2013). Simulations containing A stellar parti¬ 
cles are performed in A-body units where G = Mj = -AE^ = 1, 
in which Mt and are the total mass and energy of the clus¬ 
ter initial conditions (Heggie & Mathieu 1986). The corre¬ 
sponding dynamical time unit is 

fd = (1) 

The stellar cluster’s dynamical state relaxes over approxi¬ 
mately a relaxation time, 

^ Q.\3Wrl!^ 
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defined using the half mass radius, (Spitzer 1987; Trent! 
et al. 2010). The external tidal field of the host galaxy is im¬ 
plemented following Trent! et al. (2007b) and the reader is 
directed to section 2 of that paper for details. 

The NBODY6 distribution embeds the SSE and BSE codes 
of Hurley et al. (2000, 2002) as a simplified, on-the-fiy imple¬ 
mentation of stellar evolution^. The binaries in our clusters 
evolve with the standard BSE implementations of binary evo¬ 
lution. NBODY6 simulations include default binary evolution 
processes implemented in BSE including tidal circularization, 
magnetic braking, Roche Lobe overfiow, and binary coales¬ 
cence (Hurley et al. 2002). 

We initialize our models with an IMBH of mass Mbh at rest 
at the system center of mass. This IMBH particle then evolves 
identically to the remaining particles. In this sense, our sim¬ 
ulations do not self-consistently account for the formation of 
such an IMBH, instead we aim to simulate its subsequent dy¬ 
namical evolution in the cluster context. Because the IMBH is 
more massive than any of the stellar particles it sinks toward 
the cluster center, motivating our initial condition. We find 
within the first core relaxation time the BH acquires non-zero 
energy relative to the cluster center of mass through strong in¬ 
teractions with cluster core stars - rapidly erasing any imprint 
of this initial condition. We include gravitational wave energy 
and angular momentum losses from any very compact bina¬ 
ries following the Peters (1964) formulation, which is partic¬ 
ularly relevant for the IMBH and its companion (e.g. Miller & 
Hamilton 2002; Giiltekin et al. 2004, 2006; Leigh et al. 2014). 

Close interactions with the IMBH can lead to the tidal de¬ 
struction of a star or the inspiral and merger of a compact 
object. We handle those events as follows in our simulations. 
Eor stars, we delete any star passing within one tidal radius, 

n = (3) 

of the BH at pericenter. In a tidal disruption event, the star 
is often on a weakly-bound or near-parabolic orbit (Erank & 
Rees 1976; Lightman & Shapiro 1977). In this case approxi¬ 
mately half of the tidal debris remains bound to the BH, while 
the other half is ejected (Rees 1988). We model this behav¬ 
ior by adding half of the disrupted star’s mass to the IMBH’s 
mass following a disruption event. An ideal implementation 
of the stellar tidal disruption process would account for the 
stellar structure in determining the degree of mass loss, and 
would also allow for partial tidal disruptions (Guillochon & 
Ramirez-Ruiz 2013). Partial tidal stripping is particularly rel¬ 
evant to evolved stars which contain a dense, tidally invul¬ 
nerable core (MacLeod et al. 2012, 2013; Bogdanovic et al. 
2014). However, in this work we adopt the simplified, full- 
disruption model described above. We implement gravita¬ 
tional wave inspiral and merger events through the inclusion 
of gravitational wave corrections to the two-body equation 
of motion of tight pairs of compact objects (Peters 1964) . 
Rather than allowing merging particles to spiral to arbitrarily 
small separations, we instead merge pairs of compact objects 
whose inspiral time becomes shorter than 1 A-body time unit, 
equation (1). We do not apply a kick to the remnant Mbh, al¬ 
though such a kick may be relevant for some mergers (Holley- 
Bockelmann et al. 2008; Konstantinidis et al. 2013). 

^ Compared to the more recent stellar evolution calculations of Spera et al. 
(2015), the mass distribution of stellar remnants in our approach is favors 
lower values, especially at low metallicities, with a potential impact on the 
dynamics being the number of objects with mass much greater than the aver¬ 
age stellar mass (e.g., see Mapelli & Zampieri 2014). 


Name 

Number 

A 

Wo 

^h,0 

[pc] 

^bh,0 

[Mol 

^max 

[Gyr] 

0"^l (J ^ 

A 

3 

100k 

7 

2.3 

150 

6 

2.5 

B 

3 

100k 

7 

2.3 

150 

6 

1 

C 

3 

100k 

7 

2.3 

75 

9 

2.5 

D 

4 

200k 

7 

2.3 

150 

10.4 

2.5 


Table 1 

Table of A/^-body simulation groups A-D, the number of simulations run per 
group, as well as their parameters and initial conditions. 


We have implemented several diagnostics to allow us to bet¬ 
ter study close companions to the IMBH. We record all three- 
body encounters in which a third object approaches with peri¬ 
center distance within three times the semi-major axis of the 
most-bound companion to the IMBH. We also include the pa¬ 
rameters of all tidal disruption and merger events described 
above. Eor close tidal passages that do not lead to disruption, 
those with pericenter distances < 10^ rt we also record 

the system parameters. While we record full-cluster mass, 
position, and velocity snapshots every 10 A-body time units, 
we record the IMBH’s position and velocity every O.lA-body 
time unit. Every A-body time unit we record the parameters 
of the most tightly bound partner to the IMBH. 

2.2. Numerical Simulations 

In this paper, we present 4 groups of simulations with dif¬ 
fering initial parameters. Simulations within each group have 
statistically different realizations of the initial conditions. Ta¬ 
ble 1 defines the properties these simulation groups. We use 
the McLuster code (Kupper et al. 201 1) to generate the initial 
stellar distribution following King (1966) models with Wb = 7. 
Our simulations include a IMBH which is initialized at the 
center of mass with zero velocity and A = l-2x 10^ stars. 
The size of our clusters is constrained by the computational 
expense of the direct A-body method. These clusters are thus 
at the low mass end of the spectrum of GC masses, and a fac¬ 
tor of 10-100 less massive than typical nuclear star clusters. 
In later sections we consider how our results might scale to 
larger A clusters. 

The stars in our initial conditions follow a Kroupa (2001) 
initial mass function (IME), within the mass range O.IMq - 
3OM0, with no binaries. The metallicity is one-tenth solar. 
We initialize our clusters as tidally-underfilling by a factor of 
1.5. In A-body units, the initial tidal radius is 10.455 rather 
than 6.95 for Wb = 7 in our simulations (compare to Table 1 
of Trent! et al. 2007b). This allows for some cluster expan¬ 
sion due to mass loss before stars begin to escape. We con¬ 
sider cases of different IMBH mass, varied between 0.13% 
and 0.26% of the initial cluster mass. 

Einally, we consider two different cases for the velocity 
kicks imparted to stellar remnant NSs and BHs (for a recent 
discussion, see Contenta et al. 2015). Rather than choosing 
absolute kick velocities based on observational constraints, 
here we choose velocities relative to the cluster velocity dis¬ 
persion in order to produce different retention fractions of 
stellar remnants. This strategy is needed because the clus¬ 
ters we are able to simulate with the direct A-body method 
are small relative to realistic dense star clusters. In both cases 
remnants are given a kick drawn from a Maxwellian distri¬ 
bution with sigma of either 1 or 2.5 times the initial cluster 
velocity dispersion, These kick distributions 

imply that 19.8% or 1.6%, of kicked objects have kick veloc¬ 
ities less than which in the absence of binaries, is propor- 
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tional to the fraction of retained remnants. 

2.3. Methodological Comparison to Previous Work 

Some previous work has mentioned or studied close com¬ 
panion stars to IMBHs in GCs. Here we describe the method¬ 
ological differences in our new simulations as they compare to 
this work. Baumgardt et al. (2004b) run direct A^-body simu¬ 
lations of clusters containing IMBHs and up to 13IK stars. 
They plot several sequences of semi-major axis of tightly 
bound stars as a function of time from their simulations con¬ 
taining IMBHs (their Figures 14 and 15). These sequences 
show a long-lived most-bound companion in several cases 
with semi-major axis considerably smaller than the other stars 
and hint at the possible universality of IMBH companions. 
Baumgardt et al. (2006) focus on the production of mass¬ 
transferring binary pairs containing the IMBH and a star. 
They run A^-body simulations including approximate tidal dis¬ 
sipation and gravitational wave radiation from stellar orbits 
and study the statistics of tidally-captured stars in young (<12 
Myr) clusters as possible ULX progenitors. Konstantinidis 
et al. (2013) also model small, young (10 Myr age) clusters 
initially containing A^=32k stars using direct A/^-body integra¬ 
tion. The clusters a massive IMBH with Mbh = 500- IOOOMq . 
They study one simulation in detail, in which the IMBH cap¬ 
tures a stellar BH companion, goes through several cycles 
of exchange, and eventually inspirals and merges. Konstan¬ 
tinidis et al. (2013) argue that the merger generates sufficient 
gravitational wave recoil to eject the IMBH from their low- 
mass simulation cluster. Mapelli et al. (2013c), Mapelli & 
Zampieri (2014), and Ziosi et al. (2014) study the dynamics 
of massive stellar-mass BHs (up to 10- IOOMq) in a large 
sample of A^-body calculations of young star clusters with a 
range of metallicities and A^ = 5500 stars. 

Leigh et al. (2014) also run A^-body simulations containing 
IMBHs with between N=32k and A^=131k stars, but they do 
not include gravitational radiation from compact binaries. Af¬ 
ter the formation of an IMBH-BH compact binary, the eccen¬ 
tricity of this innermost binary varies over time. The binary 
passes many times through orbital configurations that should 
have driven it to merger. As a result, Leigh et al. (2014)’s sim¬ 
ulations show artificially long-lived companion objects, and 
a single companion persisting throughout the cluster evolu¬ 
tion. Recently, Giersz et al. (2015) have used the MOCCA 
Monte Carlo Cluster simulation code to study accretion-fed 
IMBH growth over GC lifetimes. Their simulations, particu¬ 
larly their “slow" scenario with a single massive BH, show the 
cycles of replacement of the innermost object missing from 
earlier A-body work. However, as Giersz et al. (2015) de¬ 
scribe, the Monte Carlo code is not well-suited to studying 
the detailed dynamics small number of stars tightly bound to 
the IMBH, nor is it well-suited to the inclusion of a single, 
massive object. 

Blecha et al. (2006) adopt a slightly different approach. 
Rather than simulate an entire GC, they use a binary evolution 
code to evolve the IMBH and any companions with a fixed 
background of scatterers drawn from a GC core distribution 
function. While this approach facilitates modeling a larger 
number of IMBH pairs and a detailed treatment of IMBH-star 
interactions, it does not capture the dynamical feedback of the 
IMBH and its partner on the host cluster. 

The simulations presented in this paper were designed 
specifically to study close companion objects to the IMBH in 
realistic GCs. This involved combining some of the best as¬ 
pects of the previous work described above. In particular, we 


take advantage of the accurate dynamics of the A-body meth¬ 
ods, and implement specific interaction outcomes for close 
interactions between the IMBH and stellar objects like tidal 
disruptions and gravitational wave inspirals. We model an 
IMBH which is a realistically small fraction (< 0.3%) of the 
total cluster mass. This approach limits our study to relatively 
small GCs with Aq ^ 10^, with their proportionately low-mass 
IMBHs, Mbh ^ IO^Mq, but allows us to follow their long-term 
evolution over 6-10 Gyr without approximation. In Section 4, 
we discuss some considerations in scaling the lessons learned 
from our simulations to more massive clusters or IMBHs. 

2.4. Cluster Global Evolution 

Our simulated clusters exhibit global evolution similar to 
those of others studied in detail combining IMBHs and a tidal 
field (e.g. Trent! et al. 2007a; Gill et al. 2008; Trent! et al. 
2010; Lutzgendorf et al. 2013a; Giersz et al. 2015). Thus, we 
comment only briefiy on the generic properties of our model 
clusters to set the stage for our consideration of the dynamics 
of their IMBH-hosting cores. Clusters in a tidal field slowly 
dissolve as (typically low-mass) stars escape from the sys¬ 
tem. Our models of simulation groups A and B are run for 6 
Gyr, during which time the total number of stars decreases to 

2 X 10"^. Simulation groups C and D are run to the point 
of dissolution in the tidal field, which occurs after 9 and 10.4 
Gyr, respectively. Mass loss from stellar evolution rapidly de¬ 
creases the cluster mass by ^ 30% in the first 50 Myr, then 
proceeds increasingly gradually as time goes on. The com¬ 
bination of these effects dictate that the mass range of typical 
cluster stars narrows over the course of the GC’s evolution. At 
early times, stars are the most-massive objects other than the 
IMBH, while at later times, for example > 1 Gyr, the turnoff 
mass has decreased to < 3Mq, and stellar remnants are the 
most massive objects. Lutzgendorf et al. (2013a)’s Figure 1 
gives an excellent overview of the time-evolution of A-body 
models of GCs containing IMBHs. 

The IMBH sinks nearly to the cluster center of mass, as 
for example, can be seen in Figure 2 of Konstantinidis et al. 
(2013). Once in the core, dynamical interactions with the 
IMBH and its companion object act as a heat source for the 
cluster, stalling core collapse or causing the system to expand 
(e.g. Quinlan 1996; Yu & Tremaine 2003; Trent! et al. 2007a; 
Lutzgendorf et al. 2013a). Figure 1 shows an example cluster 
evolution from simulation group A. The half-mass and core 
radii expand initially in response to stellar evolution mass loss 
and mass loss from kicked supernova remnant objects. These 
characteristic radii then stabilize for the remainder of the sim¬ 
ulation and the cluster neither expands nor contracts as en¬ 
ergy loss balances dynamical heating. Over time, the IMBH 
grows by accreting the debris of tidally disrupted stars and 
by swallowing inspiralling compact objects. The individual 
sequences are, of course, variable, but the IMBHs in our clus¬ 
ters grow by ^ 30% over the evolutionary span of the cluster 
lifetime. 

3. CLOSE COMPANIONS TO THE IMBH 

In this section, we examine the properties of close compan¬ 
ions to a cluster IMBH using our fiducial simulation group A 
described in the previous section and in Table 1 . The qualita¬ 
tive behavior of each of the simulation groups is similar, and 
we defer to Section 4 to consider differences in our results be¬ 
tween the cluster parameters represented in simulation groups 
B-D. 
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Companion 

• MS 

• Giant 

• WD 

• NS 

• BH 


time [Gyr] 

Figure 1. Characteristic radii from the IMBH in an example cluster evolution (from simulation group A). This diagram shows the evolution of the half-mass, 
Th, and core, rc, radii, along with the distance that encloses a IMBH-mass of stars, rm. Interior to these radii, we plot the semi-major axes of the 5th most-bound 
object, fl 4 , the second most bound, ai , and the most bound, gq. We plot the pericenter distance of the most-bound object, rp o, as well. Coloring of the most-bound 
object denotes its stellar-evolutionary type, labeled in the key on the right-hand side. Stellar types include main sequence (MS) and giant stars. The remnant 
types include white dwarfs (WD), neutron stars (NS) and black holes (BH). 


3.1. Stars Bound to the IMBH 

The process of capture of a close companion begins with 
the two body relaxation of the cluster. The stellar cluster’s dy¬ 
namical state relaxes over approximately a relaxation time, frh 
(Spitzer 1987). This relaxation drives the stellar distribution 
surrounding the IMBH to higher and higher densities, until a 
steady-state distribution is reached (Peebles 1972; Bahcall & 
Wolf 1976; Frank & Rees 1976; Lightman & Shapiro 1977). 
Because the IMBH mass is much greater than the average stel¬ 
lar mass, many stars can be gravitationally bound to the IMBH 
at a given time. For these stars. 
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GMbh . 1 


^bh 


+ 2''bh < 0, 


( 4 ) 


where we’ve written their orbital energy per unit mass e with 
respect to IMBH in terms of their instantaneous separation rbh 
and velocity Vbh relative to the IMBH. As such, these are the 
stars that occupy a phase space where the IMBH is energeti¬ 
cally dominant in their orbits. In general, these stars satisfy. 


a<rb = GM\yb/(jl: ( 5 ) 

where a = GMbh/26 is the semi-major axis of the bound star, 
and cTc is the core velocity dispersion. While many stars may 
exist in such bound orbits, we would not describe each of 
these as a companion to the IMBH. 

Instead, we will consider the IMBH to have a companion 
when its most-bound star is significantly more tightly bound 
than the second most bound. One way to express this is by 
comparison of those star’s semi-major axes (or equivalently 
orbital energies) with respect to the IMBH. The condition can 
be written ao/ai ^ 1 which implies so/si ^ 1. Here the co¬ 
efficient 0 denotes the most bound and coefficient i denotes 
the (i-\- l)-th most tightly bound object. In these cases, the 


most-bound object has hierarchically separated from the rest 
of the cluster and sphere of influence stars. The remaining 
stars then orbit in the combined potential of the IMBH and its 
companion. The hierarchically isolated companion can exist 
relatively unperturbed until it suffers a strong encounter. 

3.2. Companion Capture & Orbital Hardening 

How does the most-bound star end up in a orbit hierarchi¬ 
cally isolated from the other cluster stars? We explore this 
question visually in Figure 1. This diagram shows charac¬ 
teristic radii during a timeseries of one of the case A simu¬ 
lations. From top to bottom, the radii plotted are the cluster 
half-mass radius, rh, the core radius, r^, the radius enclosing 
a IMBH-mass of stars r^. Next we show the characteristic 
radii of some of the most-bound companions to the IMBH. 
These include the semi-major axis of the 5th-most bound star 
^ 4 , the 2nd-most bound star, a\, and the most bound, a^. We 
also show the pericenter distance of the most-bound object, 
^p,o- We color the most-bound object depending on its stellar¬ 
evolutionary type. 

The larger radii, those that describe the bulk of the clus¬ 
ter, like the core and the half-mass radii are nearly constant 
in time. As the radii get smaller, and describe fewer and 
fewer particles, a higher degree of variability is seen. The 
semi-major axis of the most-bound object however, exhibits 
qualitatively different behavior than the other radii, or even 
the other tightly bound objects (as shown by ai and a^). The 
most-bound star goes through a cyclic behavior of tightening 
semi-major axis followed by replacement with a new object. 
This cycle is seen to repeat throughout the timeseries of Fig¬ 
ure 1. 

The objects energetically bound to the IMBH all represent 
pairings that are ‘hard’ relative to the cluster background en¬ 
ergy (Heggie 1975; Hut & Bahcall 1983). In clusters, binaries 
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that are hard tend to get harder through interactions with re¬ 
maining cluster stars (Heggie 1975). This process can be un¬ 
derstood as a statistical trend toward energy equipartition be¬ 
tween the binary and the background as kinetic energy tends 
to transfer from the binary components to the perturber, eject¬ 
ing it faster than it entered (Hills 1975b, a). In this picture, 
each encounter leaves the IMBH plus companion pair (on av¬ 
erage) somewhat more tightly bound, and ejects the perturb¬ 
ing object with high velocity (Saslaw et al. 1974). 

Among the multiple stars bound to the IMBH, only the 
most-tightly object is ‘hard’ relative to the other bound stars. 
The orbits of less-bound objects cannot easily perturb the 
innermost orbit. The converse implies that the less-bound 
stars, despite being hard relative to the cluster core energy, are 
not hard relative to the more tightly bound companion to the 
IMBH. These less bound stars exchange energy with the in¬ 
nermost binary in a manner which prevents their orbits from 
tightening significantly. Because strong perturbations drive 
the binary hardening process, the continuous hardening of the 
most-bound star’s orbit can be approximated relative to the 
properties of the cluster background as 


da Gpc 

— (X - a 

dt (Jc 


2 

0 ? 


( 6 ) 


where pc is the mass density of the cluster core environment 
surrounding the binary (Quinlan 1996; Heggie & Hut 2003). 
This hardening rate implies that the timescale for change of 
the semi-major axis, ao/do scales as because do oc a^. 
These scalings for the time-evolution of the most-bound ob¬ 
ject’s semi-major axis suggest the qualitative behavior seen in 
Figure 1, where ao first decreases rapidly while gq is large, 
then more slowly as gq decreases. 


3.3. Companion Orbital Properties 

In this section, we examine some of the orbital properties 
that define the population of most-bound companions to the 
IMBH. Figure 2 highlights some of the key orbital demo¬ 
graphics of these companions. The lower left panel of Fig¬ 
ure 2 displays the cumulative distribution function of the log 
of ratio of the semi-major axis of the most-bound star to the 
second most bound \og{ao/a\). This diagram reveals that, in 
our simulations, the IMBH has a companion that is hierarchi¬ 
cally isolated from the rest of the cluster stars a large fraction 
of the time. 90% of the time, the companion has ao < ^ai. 
Approximately 25% of the time, ao < I0~^ai. 

The other panels of Figure 2 examine the distributions of or¬ 
bital characteristics in physical units. In the upper left panel, 
we show the orbital semi-major axis distribution of most- 
bound companions. The largest separations are around 10^ 
AU. This distance becomes similar to the average distance of 
the less-bound stars within the IMBH’s sphere of influence 
(see, for example Figure 1), and, as a result, determines the 
maximum separation at which an object is able to persist as 
most-bound without exchange. The peak of the distribution, 
lies in the range of 10^ -10^ AU, and the minimum separation 
observed in our simulations is a few AU. The upper right panel 
shows the eccentricity distribution of companion objects. It is 
nearly thermal for most of the range of e, with some decre¬ 
ment observed near ^ ^ 1, because objects that enter into or¬ 
bits that are too eccentric interact directly with the IMBH at 
orbital pericenter. The orbital separation distribution, along 
with the object masses, map to an orbital period distribution, 
shown in the lower right panel of Figure 2. 



log(ao/AU) 





logKAi) log(Po/yr) 


Figure 2. Orbital properties of closest companions to the IMBH sampled at 
equal time intervals. We plot semi major axis, ao, eccentricity, cq, and period 
Pq, of the most-bound companion. We also show the ratio, aojai, which de¬ 
scribes the compactness of the innermost companion’s orbit compared to the 
2nd most-bound star. The ao/a\ distribution shows that the IMBH has a close 
companion that is significantly more tightly bound than any other objects a 
large fraction of the time. Approximately 90% of the time, aojai < 1/3. 
25% of the time, the companion is lOOx more tightly bound than any other 
star. Typical orbital semi-major axes range from a few to a thousand AU, and 
the period distribution is also broad, peaking around 10 yr. The eccentricity 
distribution is nearly thermal, with some a slight decrement arising at ~ 1 
due to the possibility of tidal disruption by the IMBH. 


As can be inferred from an inspection of Figure 1, com¬ 
panions to the IMBH do not persist indefinitely. Instead these 
objects undergo cycles of replacement following their termi¬ 
nation - either through exchange or direct interaction with the 
IMBH. In Figure 3, we show the residence time of a given 
object as the companion to the IMBH. In the grey histogram, 
we measure the persistence time of each individual partner¬ 
ship with the IMBH. With this diagnostic, we can examine 
how long a typical close partnership between the IMBH and 
a companion lasts. Figure 3 reveals a bimodal distribution of 
persistence times of partnerships between the IMBH and other 
objects. One peak of this distribution lies around 10^-10"^ yr, 
the second peaks around 10^ yr. A tail toward longer resi¬ 
dence times connects the two peaks. The origin of the first 
peak in this distribution can be traced back directly to the or¬ 
bital period of the typical most-bound companion, as show in 
the right-hand panel of Figure 2. The longest orbital periods 
of the most-bound companion occur when ao ^ ai, and are 
10^ yr for the parameters used here. The short residence time 
branch can thus be associated with objects that persist as part¬ 
ner to the IMBH for only one orbit. Near the peak of the dis¬ 
tribution, the grey histogram in Figure 3 records objects that 
last for a small integer number of orbits. Associating the peak 
at ^ 10^ yr in Figure 3 with N orbits of of an object near the 
maximum period for most-bound companions, we find that 
^10 orbits is a typical residence time before an exchange of 
hierarchy. A tail to longer residence times represents objects 
whose orbits harden and thus persist for many more orbital 
periods. 

The origin of the secondary peak at ^ 10^-10^ yr arises 
from objects whose orbits harden sufficiently to become sig- 
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Figure 3. Residence time of closest-partner status for companions to the 
IMBH. The grey distribution shows the time span of each individual closest- 
partnership formed with the IMBH (samples per unit partnership). The blue 
distribution samples residence time of the IMBH’s parter at equal time in¬ 
tervals (samples per unit time). Thus the grey distribution shows how long 
partnerships last, while the blue distribution shows the residence time likely 
to be seen by an observer seeing the cluster at a random time. An exami¬ 
nation of the grey distribution shows that many partnerships with the IMBH 
suffer a breakup after ~ 10^ yr, a timescale associated with one to several 
orbital periods. A tail of longer-residence time objects exists, though, fol¬ 
lowed by a secondary maximum at timescales around 10^ yr, representative 
of the objects that hierarchically isolate from the remaining cluster stars. The 
blue distribution shows that these long-lasting partnerships form a dominant 
contribution to the typically observed IMBH companions. This qualitative 
conclusion can be inferred visually from Figure 1, which shows shows that 
the IMBH typically has a long-lasting companion (with residence time >10^ 
yr). 

nificantly hierarchically isolated from the remaining cluster. 
This timescale thus becomes comparable to the cluster relax¬ 
ation timescale of ^ 0.6 Gyr. In the blue histogram of Figure 
3, we show the distribution of residence times for most-bound 
companion objects sampled at even time intervals (therefore 
it is possible to count the same long-lived companion at mul¬ 
tiple sampling times in this case). If we observe the cluster 
at any given time, we are very likely to see a companion with 
residence time of 10^ -10^ yr. This impression is confirmed 
visually in Figure 1, where the populations of objects that we 
see in the timeseries last, on average, for timescales of frac¬ 
tions of a Gyr or more. 

3.4. Companion Stellar Properties 

The color coding of stellar evolutionary state in Figure 1 
shows that the demographics of the most-bound objects are re¬ 
markably diverse. This outcome is similar to the results from 
Mapelli & Zampieri (2014) for systems undergoing Roche 
Lobe overflow. Companion objects include BHs (1-2 Gyr), 
as well as MS and giant-branch stars, and WDs of various 
masses. In some cases, a single object will persist as the clos¬ 
est companion for long enough that it evolves and its stellar 
evolutionary state changes. One example of this is seen be¬ 
tween ^ 2.2-2.4 Gyr, where a single most-bound star evolves 
from the MS to the giant branch, and eventually tidally inter¬ 
acts with the IMBH. In this section, we examine the statistical 
distributions of stellar evolutionary type, mass, and radius of 
the most-bound companion objects. 

The diversity of companion types seen in Figure 1 is re- 


Name 

MS 

[%] 

Giant 

[%] 

WD 

[%] 

NS 

[%] 

BH 

[%] 

A 

40.3 

9.2 

29.9 

6.5 

14.1 

B 

11.5 

2.3 

12.5 

6.4 

67.3 

C 

27.6 

6.0 

49.6 

1.6 

15.2 

D 

21.9 

3.2 

43.6 

4.0 

27.4 


Table 2 

Demographics of stellar type of the most-bound companions to the IMBH. 
Tabulated by percentage as sampled at equal time intervals. 


fleeted quantitatively in Figure 4. The upper left panel of 
this figure shows the companions’ stellar evolutionary type. 
The most prevalent companions to the IMBH are MS stars, 
but compact objects like WDs, NSs, and BHs all represent 
significant contributions to the total demographics. The com¬ 
panion’s radius distribution, which we show in the lower left 
panel of Figure 4, directly reflects the stellar evolutionary sta¬ 
tus of the companions. We show the range of radii encom¬ 
passing WDs to giant stars in this panel. The primary peaks 
of the radius distribution lie at lO^^R©, and and > IR©, re¬ 
flecting the WD and massive MS populations, a smaller third 
peak is seen between 10^ - IO^Rq that includes stars at vari¬ 
ous phases along the giant branch. 

The objects’ mass distribution, shown in the right hand pan¬ 
els of Figure 4, exhibits a broad main peak centered at slightly 
more than a solar mass that is made up of MS stars, giants, 
massive WDs, and NSs. This peak lies substantially above the 
mean mass of cluster stars, which is initially ^ O.SSMq. We 
also observe a tail toward much higher companion masses (up 
to ^ IOMq) made up mostly of MS stars, giants, and stellar- 
mass BHs. That the mass of most of the IMBH’s companions 
is significantly larger than that of a typical cluster star is not 
surprising. The trend toward energy equipartition in two-body 
interactions causes stars more massive than the mean mass to 
sink toward orbits deep in the cluster core, where they are very 
likely to come into direct interaction (and potentially partner¬ 
ship) with the IMBH (e.g. Goswami et al. 2012; Gurkan et al. 
2004; Leigh et al. 2014). Further, more massive secondary 
objects are likely to replace an existing companion if they 
do undergo a three-body interaction (Fullerton & Hills 1982; 
Sigurdsson & Phinney 1993; Samsing et al. 2014; Mapelli & 
Zampieri 2014). In the lower right panel of Figure 4 we show 
the un-normalized distribution of companion masses in order 
to demonstrate that there is strong evolution in the compan¬ 
ion mass distribution with stellar population age. Young ages, 
< 2 Gyr, where the turnoff mass is still high, contribute the 
bulk of the companions ^ IMq . This temporal dependence 
in the companion mass can be contrasted to the approximately 
time-independent companion orbital parameter distributions 
described in the previous sub-section. 

3.5. Termination of Close Partnerships 

Close partnerships to the IMBH can be terminated through 
one of several channels. In this section, we describe the pos¬ 
sible channels and explore their relative likelihoods. 

3.5.1. Channels 

Three possible channels of termination of a close partner¬ 
ship with the IMBH exist. They are shown in the diagram of 
Figure 5 and explained below. 

• Exchange: Exchanges can occur when a third body’s 
pericenter distance from the IMBH becomes similar to 
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the semi-major axis of the binary companion gq. When 
this occurs, the system undergoes a strong three-body 
interaction, during which the gravitational attraction 
between both objects and the IMBH is similar. Either 
star may emerge from the three-body interaction as the 
new most-bound object. Statistically, there is also a 
chance that the IMBH could be the ejected object, leav¬ 
ing behind a binary of two stars. However, since the 
mass ratio is large, this exchange outcome is extremely 
unlikely (Heggie et al. 1996). Throughout the duration 
of our simulations, we record all perturbers to the cen¬ 
tral binary with pericenter distance less than 3 times gq. 
This allows us to track the statistics of these exchange 
interactions. We find that the bulk of exchanges are due 
to strong three-body interactions, where the previous 
perturber emerges from the encounter as the new com¬ 
panion. A small fraction, of order a few percent, are 
the result of more complex multi-body strong interac¬ 
tions where an object other than the strongest perturber 
emerges as the companion immediately following the 
encounter. 


Figure 4. Stellar properties of closest companions to the IMBH sampled at 
equal time intervals. The upper-left panel shows companion stellar evolution¬ 
ary type (for a key see the legend of Figure 1), the upper-right panel shows 
companion masses Mq, and the lower-left panel their radii, Rq. Companions 
include comparable fractions of stars (both MS and giants) and stellar rem¬ 
nants (WDs, NSs, and BHs). The typical mass of BH companions is substan¬ 
tially higher than the initial cluster mean mass of 0.58Mq . However, a broad 
distribution suggests that even relatively low mass companions ~ 0.3Mq oc¬ 
casionally hnd themselves tightly bound to the IMBH. The lower-right panel 
shows how the mass distribution evolves in time. At early times in the simu¬ 
lation, shown in blue, the mass distribution is weighted toward high masses, 
and includes both massive stars and massive stellar remnants. Late times, 
when the cluster turnoff mass is lower, contribute substantially to the lower- 
mass portion of the distribution. 


Companion to the 



Tidal Disruption 



time [yr] 



Channels for 
Termination of 
Partnerships 
with the MBH 




time [hr] 


Figure 5. Interaction channels that can bring about the end of a close part¬ 
nership between the IMBH and another object. Possible outcomes of a close 
partnership include three and multi-body exchanges, gravitational wave in¬ 
spirals and tidal disruptions. The bulk of partnerships are terminated through 
the exchange channel, but a fraction of long-lasting binary companions may 
end up tidally disrupted or undergoing an inspiral. If the inspiralling ob¬ 
ject is a BH or NS, the initial gravitational wave capture leads to an IMRI 
(Amaro-Seoane et al. 2007; Mandel et al. 2008; Mandel & Gair 2009; Kon- 
stantinidis et al. 2013), while if the object is a WD, tidal heating of the object 
likely leads to disruption while from an eccentric orbit and repeated flaring 
episodes (MacLeod et al. 2014). 


• Tidal Disruption: An examination of Figure 1 reveals 
that despite having the highest binding energy to the 
IMBH, the orbit of the most-bound star is still subject 
to perturbations. Not only does its orbit tighten through 
interactions with other cluster-core stars, but its orbital 
angular momentum wanders in magnitude and direc¬ 
tion, as can be seen through the orbital pericenter dis¬ 
tance, Tp 0 - As a result, the most-bound star will be oc¬ 
casionally perturbed into an orbit which leads it to pass 
within a tidal radius, rt, of the IMBH. In our simula¬ 
tions, we record these tidal disruption events and delete 
the particle from the simulation, as described in Sec¬ 
tion 2. In reality, some stars, especially those with con¬ 
densed cores, can lose part of their envelope gas with¬ 
out being completely disrupted (MacLeod et al. 2012, 
2013). The remnants of these partial tidal stripping 
events could then go on to further interaction with the 
IMBH (Hopman & Portegies Zwart 2005; Bogdanovic 
et al. 2014). 

• Gravitational Wave Inspiral: Just as the orbital wan¬ 
dering of the most-bound object permits the occasional 
tidal disruption event, if the companion to the IMBH 
is a compact object it may also undergo a intermediate 
mass ratio gravitational wave inspiral (IMRI) (Amaro- 
Seoane et al. 2007; Mandel et al. 2008; Sesana et al. 
2008; Mandel & Gair 2009; Konstantinidis et al. 2013). 
In our simulations, these are fiagged as occurring when 
the gravitational wave inspiral time becomes less than 
one A-body time unit, and the particle is removed from 
the simulation (see Section 2). Holley-Bockelmann 
et al. (2008) and Konstantinidis et al. (2013) explore 
the retention of IMBHs suffering a gravitational wave 
recoil following an IMRI event. IMBHs that merge 
with stellar mass BHs with mass ratio ^>0.1 have a 
significant probability of being ejected from the clus¬ 
ter. Smaller kicks would temporarily offset the IMBH 
from the cluster core. These inspirals typically occur 
when the companion is in a very eccentric orbit, and 
are driven primarily by the close pericenter distance, 
rather than, for example a close semi-major axis, much 
like the tidal disruptions (and in agreement with the 
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Name 

Aex 

[Gyr-‘] 

Mnsp 

[Gyr-i] 

Mde 

[Gyr-‘] 

/ex, 3b 
[%] 

/insp,0 

[%] 

/tde,0 

[%] 

A 

530.86 

1.77 

10.11 

98.62 

80.65 

43.50 

B 

253.01 

1.76 

4.54 

98.32 

83.87 

33.75 

C 

318.70 

1.76 

5.22 

99.16 

83.33 

47.89 

D 

884.06 

3.65 

11.68 

98.78 

78.95 

30.14 


Table 3 

Close encounter and termination Demographics of most-bound companions 
to the IMBH. 

findings of Gultekin et al. 2004). The A^-body sim¬ 
ulations of Baumgardt et al. (2004b) and Leigh et al. 
(2014) lack this merger channel and Leigh et al. (2014) 
demonstrate that this leads to unrealistically long-lived 
central partnerships. If the companion is a BH or NS it 
will be swallowed as it merges with the IMBH. How¬ 
ever, in some cases where the companion is a WD, 
it will tidally disrupt before merging with the IMBH 
(e.g. Casalvieri et al. 2006; Ivanov & Papaloizou 2007; 
Sesana et al. 2008; Zalamea et al. 2010; Krolik & Piran 
2011; MacLeod et al. 2014). This can proceed in one 
of two ways. In most cases, tidal heating likely leads 
to the disruption of the WD while it is still in an eccen¬ 
tric orbit. However, in some cases the orbit may fully 
circularize before the onset of Roche Lobe overflow. 

3.5.2. Termination Demographics 

The termination of partnerships with the IMBH, as well as 
other close encounters, are displayed in Figure 6. In this fig¬ 
ure, we show the same timeseries plotted in Figure 1 . Here, 
we extend the range of the plot include radii very close the 
the IMBH, including the IMBH’s Schwarzschild radius, r^. 
We also include tidal radii of typical WDs, MS stars and giant 
stars. Overset on this timeseries we highlight the orbital semi¬ 
major axis (solid symbol) and pericenter distance (open sym¬ 
bol) of encounters leading to either tidal disruption or gravita¬ 
tional wave inspiral. Tidal disruptions are shown with circular 
points colored blue for MS star disruptions and red for giant 
star disruptions. Gravitational wave inspirals are marked with 
green stars. 

Figure 6 reveals the broad demographics of changes of part¬ 
nership with the IMBH. Many changes of partnership appear 
to arise from strong interactions between the companion to 
the IMBH and other cluster stars - often including the second 
most-bound star, whose semi-major axis, a\, is occasionally 
similar to gq. When gq ^ gi changes of partnership are fre¬ 
quent, and as a result the constantly-changing companion’s 
orbit does not have the opportunity tighten significantly. An 
example of this scenario is seen in the 2.5-4.5 Gyr span 
of Figure 6. Close encounters with the IMBH often accom¬ 
pany strong encounters - many tidal disruption events either 
involve the most-bound companion or stars scattered into a 
similar orbit. The chaotic three-body dynamics of the IMBH, 
companion, and perturber system then cause a fraction of stars 
to be scattered directly toward the IMBH. An exception to this 
requirement for scatterings comes in the form of giant star dis¬ 
ruptions, which can occasionally occur when the star evolves 
to become tidally vulnerable to the IMBH rather than being 
strongly scattered (‘spoon-feeding’ the IMBH: MacLeod et al. 
2013). 

In Table 3, we examine the termination of close partner¬ 
ships with the IMBH as well as the rates of close encoun¬ 
ters with the IMBH. The most common, by far, means for 
changes of partnership with the IMBH is the exchange in¬ 


teraction, contributing approximately 98% of the changes in 
partnership. These exchanges dominate the large peak of 
IMBH companions that have residence times in the range of 
1-10^ years, as shown in Figure 3. The mean exchange rate 
is 2.5-5 X 10“^ yr“^ Of these exchanges, ^98% are sim¬ 
ple three-body exchanges in which the strongest perturber be¬ 
comes the the new most-bound companion. The remaining 
2% are more complex multi-body interactions. These may in¬ 
volve more than three objects, or they may involve the tidal 
disruption or inspiral of one or more objects. Gravitational 
wave inspirals and tidal disruption events contribute to the re¬ 
maining termination of partnerships with the IMBH. The tidal 
disruption rate is of order 10“^ yr“^ and the inspiral rate is of 
order 10“^ yr“^ in our simulations. In Section 5, we compare 
these event rates to previous predictions and mention their po¬ 
tential significance for electromagnetic and gravitational tran¬ 
sients. 

As can be inferred from Figure 6, the longer a compan¬ 
ion object persists, the more likely it is that the pair will 
be split by direct interaction with the IMBH rather than ex¬ 
change. Longer lived companion pairs segregate to smaller 
orbital semi-major axes, reducing the cross-section for en¬ 
counters or exchange in proportion to gq. Meanwhile, inter¬ 
action with the IMBH becomes progressively more likely at 
smaller gq because the orbital eccentricity needed to produce a 
close pericenter passage becomes smaller. Thus, many long- 
lived companions appear to only be replaced when they are 
scattered to their disruption, or when they undergo a gravi¬ 
tational wave inspiral. These terminations contribute to the 
large fraction of tidal disruption events and gravitational wave 
inspirals arising from the most-bound companion in Table 3. 

4. DEPENDENCE ON CLUSTER PROPERTIES 

Having outlined the basic process of IMBH companion for¬ 
mation, evolution, and destruction through reference to sim¬ 
ulation group A in the previous section, here we consider the 
dependence of our results on the chosen simulation parame¬ 
ters. We run several sets of simulations A-D in Table 1, which 
allow us to systematically explore the differences in the prop¬ 
erties of close companions to the IMBH with varying clus¬ 
ter properties. These variations prove useful in understanding 
the physics underlying companion object dynamics. Simu¬ 
lation group B has a lower kick velocity dispersion that our 
fiducial case. A, which we have up until now analyzed. This 
lower kick velocity allows a larger fraction of remnant NSs 
and BHs to be retained in or near the cluster core. Simulation 
group C has a lower IMBH mass than the fiducial case by a 
factor of two - it’s initial mass is 75Mq rather than 15OM0. 
Finally, simulation group D has a larger total number of stars 
(and mass) by a factor of 2 (2 x 10^ rather than 10^). By ex¬ 
ploring these systematic differences, we can also consider the 
extrapolation of our results to higher mass clusters or IMBHs. 

4.1. CompGfiion Objects in Cluster Evolution 

In Figure 7, we plot comparative timeseries, analogous to 
that shown in Figure 1 , for examples from simulation groups 
B-D (which may be compared directly to Figure 1 for simu¬ 
lation group A). Even looking at the examples in this figure, 
dynamical and demographic differences and similarities dis¬ 
tinguish themselves. The pattern of orbital tightening and re¬ 
placement of a most-bound companion is repeated throughout 
the panels, despite their varying cluster and IMBH parame¬ 
ters. The global evolution of the clusters is largely similar as 
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Companion 

• MS 

• Giant 

• WD 

• NS 

• BH 


Terminations 

O MS TD, Tp 

• MS TD, a 
O Giant TD, 

• Giant TD, a 
^ GW Insp, Tp 
^ GW Insp, a 
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time [Gyr] 



Figure 6. Characteristic radii measured with respect to the IMBH in the same example timeseries as Figure 1, from simulation group A. Here we expand the 
domain to include the IMBH’s Schwarzschild radius, rs, and estimated tidal radii of typical WDs (brown), MS stars (blue), and giants (red). Overset on this 
timeseries we plot the semi-major axis (solid symbols) and pericenter distance (open symbols) of tidal disruption events (circular points) and gravitational wave 
inspirals (stars). These direct interactions with the IMBH often lead to the termination of long-standing partnerships with the IMBH. On visual inspection, several 
categories of events emerge. Gravitational inspirals are particularly likely to come at the end of a long partnership when the orbital SMA has decreased. Tidal 
disruptions, by contrast, often come in times of multi-body interaction that scatter stars toward the IMBH. These are particularly likely to occur when ao ^ ai, 
either because ao is large or is smaller than average. In several cases represented here a most-bound star evolves off the MS and up the giant branch until it 
eventually reaches the point of tidal disruption due to its growing radius. 


well: cluster core and half mass radii are relatively constant 
in time, just as with simulation group A. 

Differences in the stellar dynamics and demographics of the 
companion objects also become apparent on closer inspection 
of Figure 7. Simulation group B, with a larger population 
of retained BHs, shows a larger quantity of BHs joining the 
IMBH as most-bound companions. The cycles of orbital hard¬ 
ening experienced by these BHs are longer-lived than those 
of most of their stellar counterparts. By contrast, simulation 
group D, with a larger total number of stars than the fiducial 
case shows more frequent exchanges of partnership, but a di¬ 
versity of companion-object stellar evolutionary types. 

These BHs preferentially enter into partnerships with the 
IMBH because of their significantly higher-than-average 
mass, which causes them to sink in the cluster potential due 
to dynamical friction. When these stellar-remnant BHs inter¬ 
act with the IMBH and its companion, their high mass makes 
them very likely to replace an existing companion. Heggie 
et al. (1996) has shown that the exchange cross section de¬ 
pends quite strongly on the mass of the tertiary, perturbing 
body, even in the case where the mass ratio between the binary 
components is large. In Figure 8, we plot the fraction of en¬ 
counters that result in an exchange given a particular perturber 
to companion mass ratio, M^/Mq. This fraction is computed 
by counting the fraction of encounters in which the perturber 
comes within ao of the binary center of mass at pericenter that 
result in an exchange. While the exchange fraction depends 
strongly on perturber mass ratio, it does not depend nearly 
as strongly on other cluster properties. This plot implies that 
massive objects can exchange into the inner binary with rel¬ 
ative ease and that, as companions, they can undergo many 


encounters, each hardening their orbit before one leads to an 
exchange. 

4.2. Companion Demographics 

Figure 9 further traces how companion demographics prop¬ 
agate to orbital distributions. The histogram of orbital semi¬ 
major axis is shown for each of the simulation groups A-D, 
and is broken down by companion type. Visually from Figure 
7, and here quantitatively, we can infer that types of com¬ 
panions have different sub-distributions. When combined in 
varying percentages, these demographic differences lead to 
variations in the overall semi-major axis distribution of IMBH 
companions. The most obvious example are BHs. Based on 
the histograms of Figure 9, BHs exist in orbits which are pref¬ 
erentially among the most tightly bound. The distribution is 
broad and over the course of the simulation there were BHs 
at all separations. We can understand this distinction because 
BHs are preferentially difficult to replace in exchange interac¬ 
tions due to their high mass, and are not prone to tidal disrup¬ 
tion when the the companion has orbit with e ^ 

By contrast, giant stars tend to exist in the least-bound or¬ 
bits in Figure 9. These objects are especially tidally vul¬ 
nerable, and so are subject to tidal disruption by the IMBH 
at larger orbital separations for a given eccentricity than are 
more compact MS stars or WDs. The cases of MS stars and 
WDs lie intermediate to that of giants and BHs. MS stars’ in¬ 
nermost semi-major axes are also limited by tidal disruption 
as their orbits cycle through high-eccentricity phases (note the 
rapidly varying eccentricity of the innermost orbit in Figure 
1). WDs are less tidally vulnerable, but are not as massive as 
the BHs, so their distribution peaks at somewhat larger sepa- 
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Figure 7. Comparative timeseries of first 6 Gyr of cluster evolution of examples from simulation groups B-D. This figure is analogous to Figure 1, which plots 
the evolution of one of the simulation group A models. As in Figure 1, we color hard-binary companions to the IMBH based on their evolutionary type. Table 1 
describes the differences between simulation groups B-D as compared to A. 


rations. The case of the NSs is interesting and somewhat more 
subtle. Although these remnants are ^ 3 x the average stellar 
mass in the cluster, they are subject to kicks on formation and 
they do not re-segregate to the cluster core as efficiently as the 
more-massive BHs. 

Figure 10 quantitatively explores how these demographic 
differences propagate to the other properties of the most- 
bound companions to the IMBH. The upper left panel shows 
companion object demographics, which are also tabulated in 
Table 2. Qualitatively, the demographics remain unchanged in 
that the companions include stellar remnants like WDs, BHs, 
and NSs, along with MS and giant stars. As seen in Figures 
7 and 9, the most dramatic difference arises with simulation 
group B, and the retention of a higher fraction of BHs in the 
cluster’s central regions. Simulation group D, with its larger 
total particle number, also exhibits a mild enhancement in the 


fraction of BH companions. This difference can be traced 
to the larger total number of BHs retained in these simula¬ 
tions (despite the small retention fraction). Companion de¬ 
mographics shape the companion mass distribution, shown in 
the upper right panel of Figure 10. The enhanced BH fraction 
is refiected in the higher peak at companion masses > SMq 
for simulation groups B and D. 

4.3. Companion Dynamics 

In addition to demographic changes in the typical proper¬ 
ties of IMBH companions, there are also dynamical differ¬ 
ences we can observe across our simulations. We begin by 
examining how the long-lasting tail of the residence time dis¬ 
tribution of IMBH companions is affected by changing cluster 
properties in the central panel of Figure 10 (the residence time 
plotted is counted per-partnership as is the grey distribution in 
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log(Mp/Mo) 

Figure 8. Fraction of close encounters between perturbers of mass Mp and 
the IMBH and companion pair that result in an exchange. The exchange 
fraction depends sensitively on the perturber to companion mass ratio, al¬ 
lowing massive objects to exchange into the inner binary easily, and then to 
resist being exchanged out by subsequent encounters with lighter objects. By 
contrast, we see little dependence of the exchange fraction on other cluster 
properties, like the IMBH mass. Errorbars are poisson errors on individual 
bins in the exchange cross section. 



0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 


log(ao/AU) log(ao/AU) 

Figure 9. Distributions of companion object semi-major axis, ao, subdivided 
to show companion demographics. Panels show simulation groups A-D out¬ 
lined in Table 1. These panels illustrate how differing companion demo¬ 
graphics between simulation groups propagate to different distributions of 
companion orbital properties. The distributions of compact object (BH and 
WD especially) semi-major axes peak at smaller values than the overall dis¬ 
tribution, so varying quantities of compact object companions propagate to 
different overall distributions of companion orbits. 


Figure 3). This tail of long-lived companion objects is shaped 
by processes of companion exchange and disruption. As a re¬ 
sult, it is influenced by a combination of cluster, IMBH, and 
companion demographic properties. An increased fraction of 
BH companions (as in simulation case B) implies more long- 
lasting partnerships because the IMBH’s companion is con¬ 
siderably more immune to exchange interactions (see Figure 
8) and to tidal disruption. In fact, we can observe this directly 
in the termination statistics of Table 3, where we see that ex¬ 
changes are significantly less frequent in simulation group B 



log(tres/yr) 



Figure 10. In this figure, we show a general comparison of IMBH compan¬ 
ion properties with varying simulation initial parameters. We include four 
groups of simulations in this diagram, A-D, described in Table 1. Simulation 
group B, has a lower kick velocity than group A, as a result, more stellar 
remnant BHs remain in the cluster core. These BHs are more massive than 
the typical stellar companion, and therefore have longer residence times in 
companionship with the IMBH because more massive innermost binaries are 
more difficult to disrupt. As a result of the increased time for the binary 
to harden, the ao/ai distribution extends to smaller separations. Simulation 
group C includes a lower IMBH mass of 15Mq instead of I50Mq, this has 
only a mild effect on the stellar properties of companions, but does allow 
companions to segregate to tighter separations in and thus also shorter 

median orbital periods. Finally, simulation group D includes twice as many 
initial particles as the other simulations, N = 200k. The increased particle 
number implies shorter residence times as the IMBH’s most-bound compan¬ 
ion because there are more potential perturbing bodies, and also shorter or¬ 
bital periods because there are, on average, a larger number of tightly-bound 
stars. 


(because of the mass sensitivity of the exchange cross-section. 
Figure 8). Tidal disruptions are also nearly a factor of two less 
common than in case A, with a smaller fraction arising from 
the most-bound companion. With lower IMBH mass, as eval¬ 
uated in simulation group C, weaker gravitational focusing of 
surrounding stars leads to the lower exchange rate tabulated 
in Table 3. This propagates to more objects with long res¬ 
idence times in the 10^ - 10^ yr range. The larger number 
of stars present in simulation group D, by contrast, implies a 
larger encounter and exchange rate and shorter typical resi¬ 
dence times. 

The lower panels of Figure 10 explore differences in com¬ 
panion periods, Pq, and hierarchies, aojai. Here distinctions 
arise both in the form of companion demographics and IMBH 
mass. All of simulation groups B-D exhibit shorter typical 
orbital periods than simulation group A. The peak in the or¬ 
bital period distribution shifts by about 1 dex from ^ 10 yr, to 
^ 1 yr, but the distribution remains broad in all cases. Com¬ 
panion hierarchy, as measured by ao/^^i, is also affected both 
by companion demographics and IMBH mass. The flducial 
case, simulation group A, shows companions which are less 
segregated than groups B-D. In the case of simulation group 
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B, the tighter companion distribution can be traced to the 
increased BH fraction and correspondingly longer compan¬ 
ion residence time (which allows the IMBH-companion pair 
more time to harden prior to disruption). Simulation group 
D’s tightly-bound tail can be similarly explained by a mildly- 
increased fraction of BH companions. 

The IMBH mass effect on companion hierarchy can also 
be observed in the lower-left panel of Figure 10. The median 
ai^jai shifts by a factor of approximately 3 between simula¬ 
tion groups A and C, which have a factor of 2 difference in 
IMBH mass. The origin of this distinction can be understood 
through a consideration of the dynamics of orbital hardening 
and exchange for IMBH companions. A typical close pas¬ 
sage of a perturbing body of mass Mp through the IMBH- 
companion pair carries away an energy of order 


A£p 


GMpMp 

2^0 


(7) 


when it passes with orbital pericenter within the semi-major 
axis of the inner binary (Yu & Tremaine 2003). As com¬ 
pared to the orbital energy, £"0 = GMbhMo/2(2o, this typical 
per-encounter change of energy is 


A£p _ ^ 
^0 ^bh 


( 8 ) 


So the fractional, per-encounter hardening rate of the IMBH 
and companion pair scales inversely with the IMBH mass. 
But, the rate of encounters also changes with IMBH mass. 
The encounter cross section is dictated by gravitational focus¬ 
ing onto the inner binary, and is Ygnc ~ 27rGMbh^2o/v^- Asso¬ 
ciating a typical velocity with the cluster core velocity disper¬ 
sion, the encounter rate, Aenc where is the core 

number density, thus scales as. 


Aenc OC^cCTc^^bh^^O, 


(9) 


proportional to both Mbh and ap, as long as the low-angular 
momentum phase space of encountering orbits is efficiently 
replenished (e.g. Frank & Rees 1976). As a result, although 
each encounter with a perturber carries away less energy 
oc , the encounters are more common oc Mbh and the inner¬ 
most binary hardens at an identical rate regardless of IMBH 
mass. 

Although the hardening rate remains unchanged with dif¬ 
fering IMBH masses, the residence time, as shown in the cen¬ 
tral panel of Figure 10, does change. In particular, when the 
encounter rate goes up, so does the exchange rate of IMBH 
companions. This is because the exchange fraction in our sim¬ 
ulations does not appear to depend strongly on IMBH mass, 
as shown in Figure 8, instead it is sensitive primarily to the 
perturber-companion mass ratio. As IMBH masses go up, 
their companions harden at the same rate, but are exchanged 
with increasing frequency, restarting a new hardening cy¬ 
cle. This result implies, when extrapolated to larger IMBH 
masses, that only IMBHs of moderate mass ratio to their sur¬ 
rounding stars (perhaps those in the range r\j 
will ever exhibit tightly segregated companions. 


5. DISCUSSION 

In this section we draw on our results concerning the prop¬ 
erties of IMBH companions to explore how these objects 
can occasionally reveal the presence of an otherwise invisi¬ 
ble IMBH in the cluster core. 


5.1. IMBH Companions in Old Clusters 

We have focused so far on IMBH companions across a 
range of stellar system ages. However, most GCs are very old 
stellar systems, with ages >10 Gyr (e.g. Marm-Franch et al. 
2009; Forbes et al. 2015; Trent! et al. 2015). How might the 
companions to IMBHs in these very old stellar systems dif¬ 
fer? To examine this question we look at several evolutionary 
sequences of the last several Gyr of cluster evolution of the 
models in simulation group D. These evolutionary tracks are 
shown in Figure 11 . We find that by these late epochs inter¬ 
actions with the IMBH (either mergers or ejections) have ex¬ 
hausted the supply of stellar-mass BHs which might fall into 
partnership with the IMBH (Leigh et al. 2014). If realistic 
GCs begin with much larger masses and A than their present- 
day properties (or our smaller-still simulated systems), they 
would retain a similar fraction but larger number of BHs to 
late times (Morscher et al. 2015). If these BHs remain in 
the cluster we would expect them to mass-segregate toward 
the cluster core and, perhaps, play a role in shaping the time- 
dependence of the companion mass function, show in Figure 
4. At 10 Gyr, the turnoff mass is now ^ IMq, and is there¬ 
fore less than the mass of typical stellar remnants like WDs 
and NSs. As a result, the most common companions for the 
IMBH in our simulations at these times are massive WDs, 
those with masses > IMq , which are therefore more massive 
than the MS stars in the cluster core. NSs similarly form a 
fraction of late-time companions, because those that were re¬ 
tained in the cluster (but perhaps kicked to the outskirts) have 
had sufficient time to re-segregate into the cluster core. 

5.2. Dynamical Mechanisms for Repeated Tidal Disruption 
Flares by IMBHs 

One motivation for performing this study was to access 
whether close companions to the IMBH are ever driven into 
multiple-passage tidal interactions. Such multi-passage tidal 
disruption fiaring interactions can be difficult to achieve. For 
multiple passages of similar strength to be sustained, the star 
must remain bound to the IMBH with its orbit relatively un¬ 
changed (in energy, or more importantly, angular momen¬ 
tum). Further, the object undergoing the passages must ex¬ 
pand upon mass loss either through its adiabatic response to 
mass-loss or heating due to the violent interaction. These re¬ 
quirements suggest that orbits leading to multi-passage dis¬ 
ruptions must exist in a phase space where the per-orbit scat¬ 
ter due to the remaining cluster stars is small (MacLeod et al. 
2013). For MS stars or WDs in clusters with as few stars as we 
consider here, the only star in such an orbit is the most-bound 
companion to the IMBH. Thus, if this most-bound star enters 
into a weak tidal encounter with the IMBH the stripping may 
persist for several passages. 

In our simulations, we find that a reasonable fraction (30- 
50%) of tidal disruption events involve the most-bound star 
(Table 3). Although we track only full disruptions dynami¬ 
cally in our simulations, many of these would, in fact lead to 
partial tidal disruptions over multiple orbits. Stars with con¬ 
densed cores relative to their envelopes, in particular evolved 
stars and those with radiative envelopes are particularly able 
to retain a bound core despite losing part of their envelope 
(MacLeod et al. 2012; Guillochon & Ramirez-Ruiz 2013). 
We find that some giant stars evolve to the point of mass trans¬ 
fer with the IMBH (MacLeod et al. 2013). WDs that undergo 
gravitational wave inspiral are also very likely to disrupt be¬ 
fore their orbits circularize, leading to episodic mass transfer 
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time [Gyr] 


Figure 11. Evolution of representative members the four simulation group D models from 6 Gyr to their tidal dissolution at ~ 10.5 Gyr. The lines plotted here 
are a simplified set of those presented in Figure 1, from bottom to top in each panel: ao (color), ai (grey), rm, ^c, rh. This Figure may be compared with the first 
6 Gyr of evolution shown in Figures 1 and 7. IMBH companions in old clusters are qualitatively similar to those at earlier epochs, but contain fewer BHs and a 
larger fraction of WDs. Sequences for other realizations of a given simulation group are quantitatively different but qualitatively similar to the models reproduced 
here. 


and repeated flaring (MacLeod et al. 2014). Many of the dis¬ 
rupted most-bound stars reach the tidal radius through a secu¬ 
lar interaction with other tightly-bound objects - opening the 
possibility of a large number of passages with similar pericen- 
ter distance. In the ensuing repeated tidal stripping episodes 
the individual flares would be qualitatively similar to other 
tidal disruption flares albeit with less total fluence. The rep¬ 
etition time between flares would be governed by the orbital 
period distribution of the most-bound companion. 

This leads us to the case of the extraordinary repeated- 
flaring transient HLX-1 (Farrell et al. 2009) by a putative 
IMBH with Mbh ^ IO^Mq (e.g. Davis et al. 2011). If this 
object is powered by the partial tidal stripping of a star in an 
eccentric orbit the object must have gone through >10 pas¬ 
sages based on the lightcurve, with orbital period near a year. 
The IMBH mass predicted for HLX-1 is nearly a factor of 100 
larger than those we are able to simulate. However, with the 
scalings of section 4 as guidance, we can speculate whether 


these properties appear plausible given our results. 

First, having tightly bound stars in orbit around the IMBH 
appears likely given the IMBH (rather than supermassive) na¬ 
ture of the putative HLX-1 BH. Although scalings to a more 
massive system would modify them, our distributions of or¬ 
bital periods easily extend to a orbital periods in the years to 
tens of years range. From a pericenter distance slightly out¬ 
side the tidal radius, the donor star could have been tidally 
excited until it started transferring mass to the IMBH (Hop- 
man et al. 2004; Baumgardt et al. 2006; Blecha et al. 2006). 
Secondly, a large fraction of tidal disruption events in our low- 
density, small-A clusters originate from the most-bound com¬ 
panion to the IMBH (as compared to the supermassive BH in 
a galactic center context upon which most tidal disruption lit¬ 
erature is deflned). We observe that the most-bound stars that 
tend to interact with the IMBH in our simulations arrive in 
disruptive orbits through secular interaction, rather than sin¬ 
gle, strong scattering events. Thus, it would not be entirely 
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surprising if a tightly-bound star in the HLX-1 system ended 
up in a grazing, rather than fully-disruptive, interaction with 
the IMBH. In conclusion, a scenario involving a hierarchi¬ 
cally isolated companion star being driven to interact with the 
IMBH is a scenario that qualitatively explains the properties 
of a repeated flaring source like HLX-1, and extrapolation of 
our results indicates that such a conflguration can be estab¬ 
lished through dynamical interactions of stars with the IMBH. 

5 . 3 . Revealing IMBHs in GCs 
5 . 3 . 1 . Luminous Companions to a Dark IMBH 

The calculations of this paper have demonstrated that 
IMBHs residing in clusters acquire and retain close compan¬ 
ions. In rare periods in which the IMBH disrupts and accretes 
its companion, a luminous accretion signature will mark the 
presence of the IMBH. We discuss this possibility further in 
the following section, but here we consider the ability of com¬ 
panion objects to reveal the IMBH in periods of relative qui¬ 
escence. Among the companion objects, approximately half 
are luminous stars like MS stars and giants, while the other 
half are WDs, NSs and smaller BHs. Can this luminous pop¬ 
ulation of stars serve to reveal the presence of otherwise dark 
IMBHs in nearby clusters? 

First, the binary nature of the IMBH plus companion pair 
could potentially reveal the IMBH’s presence. The typical 
companion orbits found in our simulations (semi-major axes 
of ^ 10^ AU and periods of ^ 10 yr, shown in Figure 2) imply 
typical orbital velocities of hundreds of km s“^. If spectro¬ 
scopically detected, this orbital motion would strongly point 
to a massive companion because of the high velocity given 
combined with relatively long orbital period (compared to a 
stellar-mass binary). This approach extends work already be¬ 
ing done to explore radial velocities of stars bound to potential 
IMBHs in cluster cores (e.g. Noyola et al. 2010; Lutzgendorf 
et al. 2013b). But we argue that among the tightly bound stars, 
it would be reasonable to expect a single very compact com¬ 
panion orbit in a large fraction of clusters hosting a IMBH. 

Luminous stars also shed substantial stellar winds. With 
the star tightly bound to the IMBH we examine whether 
and how these winds might contribute to fueling accretion 
activity onto the IMBH. The relevant scale for compari¬ 
son here is the Bondi accretion rate at which the IMBH 
could accrete gas from the cluster-core interstellar medium, 
^bondi = 47r(GMbh)^p/Cs (Volonteri et al. 2011). This inter¬ 
stellar medium is shaped by the interacting stellar winds of 
all the stars in the cluster core. Taking typical values of 
10“^^ g cm“^ and sound speed = 10^-10^ cm s“\ we And 
4^bondi ~ 1O“^^-1O“^^M0 yr“^ (e.g. Naimanetal. 2012,2013). 

In Figure 12, we compare this nominal accretion rate from 
the background of stars to that fed by the companion. Here 
we plot the duty cycle (fraction of time) spent accreting 
above a given level. The upper series of lines in this flg- 
ure shows Mq, the mass loss rate of the companion star, es¬ 
timated using Reimers (1975) mass loss formula, M = 4 x 
lO“^^7^/?(LoMo/Ro)/(^04^q/^0) with pr = 0.5 (McDonald & 
Zijlstra 2015). We then estimate the fraction of this mass lost 
which is bound to the IMBH using the formulae of Miller 
et al. (2014, described in their section 4.1) The bound fraction 
depends on the orbital properties through the semi-major axis 
and eccentricity, which determine the velocity of the wind- 
shedding star relative to the IMBH. It also depends on the 
terminal velocity of the winds themselves, which we take to 



Figure 12. The distribution of accretion rates fed by the stellar wind of the 
companion to the IMBH. The series of lines labeled Mq is the wind mass loss 
rate of IMBH companions estimated with Reimers (I975)’s formula. The 
lower series of lines use the wind velocity and orbital properties to estimate 
the fraction captured by the IMBH. The shaded region denotes the range of 
possible Bondi accretion rates onto a 150 Mq IMBH with ambient gas density 
of 10“^^ g cm“^ and sound speed Cs = 10^ -10^ cm s“^ (Naiman et al. 2012, 
2013), where Mbondi = 47r(GMbh)^p/cs. Stellar winds are only captured by 
the IMBH less than 10% of the time. The typical wind-fed accretion rates 
during this 5% duty cycle are competitive with Bondi accretion from the 
background of cluster-core stars but only strongly dominant if our more con¬ 
servative estimate of the gas properties (cs ~ 10^ cm s“^ ) applies. Thus, we 
do not expect companions to IMBHs to impose strong limits on the accretion- 
fed luminosity of cluster IMBHs. 

be the star’s escape speed, ^JiGMq/Rq. The slow, massive 
winds of giant stars, for example, are particularly easy to cap¬ 
ture. We plot the portion of companion wind material that 
is bound to the IMBH as Mbh in Figure 12. Most of these 
captured winds come when the IMBH hosts a giant star com¬ 
panion. These cumulative distributions show that the IMBH 
captures the winds of its companion at rates comparable to 
4^bondi of Order 5% of the time in our clusters. 

As a result, we can expect that the > 50% of stellar IMBH 
companions should illuminate their host BH through their 
winds, providing a minimum accretion rate onto the IMBH. 
However, the fact that the IMBH has a companion does not 
imply that most IMBHs should accrete substantially above the 
Bondi accretion rate from the ambient medium. Of course, 
how these alternative fuel sources transfer into accretion lu¬ 
minosity represents a signiflcant uncertainty. Media captured 
from the cluster core has low net angular momentum, and 
likely does not efficiently form a disk (e.g. Cuadra et al. 2008), 
perhaps contributing to a low radiative efficiency. There is 
still more uncertainty in how much of the wind-deposited 
gas might reach a IMBH in scenarios with realistic feedback 
(Naiman et al. 2009; Park & Ricotti 2011, 2012, 2013) and 
what fraction of the resultant accretion energy would con¬ 
tribute to radio and X-ray luminosities (Maccarone 2004; 
Maccarone et al. 2005; Strader et al. 2012; Wrobel et al. 
2015). Given these uncertainties, wind feeding from the com¬ 
panion implies a minimum activity level for the fraction of 
IMBHs hosting giant star or MS companions. 

5 . 3 . 2 . Electromagnetic and Gravitational Transients 
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Stellar tidal disruption events lead to a stream of debris 
from the disrupted star accreting onto the IMBH. The ensu¬ 
ing accretion flare may be extremely luminous (Ramirez-Ruiz 
& Rosswog 2009). For the relatively low-mass clusters with 
low-mass IMBHs we have considered through our direct N- 
body integrations, the rate of tidal disruptions is relatively 
low, ^ 10“^ yr“^. More massive IMBHs and denser clus¬ 
ters both imply higher event rates, with A^tde oc (Rees 

1988). The events themselves, however, are likely short, with 
the flare characteristic timescales of order weeks to months 
(e.g. Ramirez-Ruiz & Rosswog 2009; MacLeod et al. 2012; 
Guillochon & Ramirez-Ruiz 2013) and observable signatures 
lasting at most for years or tens of years through nebular 
emission (Clausen & Eracleous 2011; Clausen et al. 2012). 
The implied duty cycle of these flaring episodes is thus ex¬ 
tremely small and the fraction of cluster IMBHs expected to 
be in a flaring state is ^ 10“^. Tracing this number of clus¬ 
ters is not inconceivable, optical surveys for tidal disruption 
flares routinely cover ^ 10^ galaxies (e.g. van Velzen et al. 
2011), each of which might host 10^-10^ globular clusters 
(Brodie & Strader 2006). But this scaling argument is ham¬ 
pered by the fact that the Eddington-limited luminosity of low 
mass IMBHs is much lower than their supermassive counter¬ 
parts. Thus, this argument only suggests a reasonable fraction 
of detectable events if the occupation fraction of IMBHs in 
local-universe clusters is of order unity. 

Relativistically beamed jet emission from tidal disruption 
events may be visible in some cases, and for low mass BHs, 
the implied luminosities are large enough to suggest that these 
events could be observed to cosmological distances (Bloom 
et al. 2011; Zauderer et al. 2011; Berger et al. 2012; De Colle 
et al. 2012; Zauderer et al. 2013). Because of their high appar¬ 
ent luminosity these high-energy transients, especially those 
arising from rare but luminous WD tidal disruptions, could 
dominate the observable population of tidal disruption flaring 
events arising from IMBHs (Krolik & Piran 2011; MacLeod 
et al. 2014, 2015). With a population of these transients we 
might eventually hope to understand the cosmological demo¬ 
graphics and temporal evolution of IMBHs in clusters. But 
given their extreme rarity, these events do little to probe any 
potential local-group population of IMBHs in CCs. 

Cravitational wave inspirals do offer a promising new probe 
of the presence or absence of IMBHs in nearby clusters (e.g. 
Miller & Hamilton 2002; Cultekin et al. 2004; Hopman & 
Portegies Zwart 2005; Amaro-Seoane et al. 2007; Ivanov & 
Papaloizou 2007; Sesana et al. 2008; Mandel et al. 2008; 
Mandel & Gair 2009; Konstantinidis et al. 2013). With laser 
interferometer experiments like advanced LIGO entering a 
phase of enhanced sensitivity, the possible detection of grav¬ 
itational radiation from a merging binary appears promising. 
How might gravitational wave inspirals of compact objects 
into cluster IMBHs compare? Because of their higher masses, 
IMRI events enter the LIGO band only near the end of the 
inspiral (Miller & Hamilton 2002). For example, the grav¬ 
itational wave frequency at the innermost stable circular or¬ 
bit of a IOOMq IMBH is ~ 40 Hz, with an amplitude that 
might allow its observation to z ^ 0.1 (Mandel et al. 2008). 
Similarly, Konstantinidis et al. (2013) show that a LISA-like 
space-based detector could capture lower-frequency emission 
from IMRIs at times several years prior merger out to z ~ 0.7. 

The event rate of IMRIs compared to standard binary merg¬ 
ers is much more uncertain. In our simulations, we And simi¬ 
lar IMRI rates per GC, ^ 1 Gyr“\ as estimated previously by 


Mandel et al. (2008). For context, we can assume a galaxy 
density of 10^ Gpc“^ and 10^ globular clusters per galaxy (Si- 
jacki et al. 2015; Brodie & Strader 2006). If we conservatively 
take the inspiral rate from our low-mass clusters as represen¬ 
tative, then an inspiral rate of order 1 yr“^ Gpc“^ is expected 
if every cluster hosts a IMBH. This is a factor of ^ 100 lower 
than the double NS inspiral rate estimated recently (e.g. Do- 
minik et al. 2013, 2015), but the large chirp masses of these 
IMBH binaries make their inspirals observable in the larger 
volumes mentioned above that might imply similar detection 
rates of both classes of source (Miller & Hamilton 2002; Man- 
del et al. 2008) These constraints are compelling enough to 
offer hope that the upcoming era of gravitational wave de¬ 
tections will either offer guiding detections or constraining 
non-detections to pin down the IMBH occupation fraction in 
low-redshift clusters. 


5.4. Extensions and Future Work 


We have explored the question of close companion ob¬ 
jects to IMBHs in dense cluster environments using direct 
A-body numerical simulations in this paper. This chosen ap¬ 
proach carries the advantage of direct integration of each par¬ 
ticle’s equation of motion, ensuring accurate dynamics. How¬ 
ever, it necessitates some approximations due to the compu¬ 
tational expense of each individual simulation. In particu¬ 
lar, in our analysis we consider clusters with initial mass of 
^ 5 X IO^Mq or IO^Mq, by the time these clusters evolve, 
they are less massive (and thus less dense) than the the typi¬ 
cal Milky Way GC. We similarly include IMBHs of relatively 
low initial mass, 75 or 15OM0. More massive clusters might 
harbor more-massive still IMBHs, changing the IMBH to star 
mass ratio. While these constraints are shared by every di¬ 
rect A-body simulation method, they necessitate a consider¬ 
ation of how the stellar dynamics will scale and extrapolate 
to more massive systems found in Nature. We have qualita¬ 
tively explored the response and dependence of our results to 
dependencies in cluster mass and density and IMBH mass in 
Section 4. To extend our simulated dynamics into the regimes 
that found in more massive IMBH environments, an A-body 
algorithm particularly well suited to study more extreme mass 
ratios, like that recently published by Karl et al. (2015), may 
offer a promising numerical approach. 

We have included no primordial binaries in the simulations 
presented in this work. This simpliflcation is not realistic; 
binary fractions are thought to be of order a few percent in 
typical GCs (Heggie & Hut 2003). In future work, we intend 
to include a realistic binary fraction and study the role these 
primordial binaries may play in interacting with the IMBH 
and its close companions (Pfahl 2005). Before addressing 
this question with numerical simulations, we can make a few 
analytic estimates of their potential importance. In a cluster 
environment, only hard binaries survive the interactive clus¬ 
ter core environment (Heggie 1975). To order of magnitude, 
these binaries have separation GMbm/cr^. If such a binary 
passes close to the IMBH and is split by the Hills (1988) 
mechanism, the one star can remain bound to the IMBH while 
the other is ejected. The bound star has a minimum typical 
semi-major axis of 


^min ^ 



2/3 

^bin? 


( 10 ) 


plugging in a few typical numbers for our simulation clus¬ 
ters, this might indicate ^ 3 x 10^^ cm, for marginally- 
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hard source binaries. This is near the peak of the semi-major 
axis distribution of IMBHs companions in our case A sim¬ 
ulations (as seen, for example, in Figure 2). Thus, in some 
cases, these split binaries would leave behind stars with sim¬ 
ilar semi-major axis to the most-bound IMBH companion. 
Additional strong interactions between the newly-bound star 
and the companion would ensue with the possibility of an ex¬ 
change of hierarchy. 

Trenti et al. (2007a) run direct N-body simulations of star 
clusters with primordial binaries and a central IMBH under 
idealized conditions (N = 16384 equal-mass particles) and 
find that the IMBH presence increases the binary disruption 
rate in qualitative agreement with the theoretical predictions 
by Pfahl (2005). However, the large majority of the binaries 
are split by interactions with other stars in the dense central 
stellar cusp surrounding the IMBH, rather than by the IMBH 
tidal field (Trenti et al. 2007a, Figure 6). This hints that a 
majority of binaries in relatively low angular momentum or¬ 
bits may be split before reaching the IMBH, and that those 
that survive as bound pairs and experience a close encounter 
with the IMBH are likely to be preferentially harder binaries. 
The earlier idealized experiments of IMBH-binary interac¬ 
tions also highlight the importance of resolving the dynamics 
of the sphere of infiuence exactly through direct N-body sim¬ 
ulations. We plan to explore these effects more fully with the 
aid of future simulations. 

6. SUMMARY & CONCLUSION 

IMBHs residing in clusters acquire close companions 
through strong dynamical interactions with their surround¬ 
ings. We use full A-body calculations including stellar evo¬ 
lution to present a systematic study of the demographics, dy¬ 
namics, and observable properties of these IMBH compan¬ 
ions. This work confirms suggestions that IMBHs in dense 
clusters should have close companions for the majority of 
their evolution (Blecha et al. 2006; Mapelli et al. 2013c; 
Mapelli & Zampieri 2014). We find that ^ 90% of the time, 
companion objects have semi major axes less than a third that 
of the next-most bound object, ao <1 /3ai . The hierarchy can 
be substantially more extreme, with a few percent of compan¬ 
ion configurations having ao (Figure 2). The most- 

bound object demographics are broad, including a wide vari¬ 
ety of stars and compact objects. 

We show that most-bound partners to the IMBH go through 
cycles of orbital hardening followed by destruction, illustrated 
in Figure 1. The typically-observed pairings are long-lived, 
with residence times of 10^ -10^ yr, shown in Figure 3. These 
pairings are terminated through exchange or close interaction 
with the IMBH (gravitational wave inspiral or tidal disrup¬ 
tion). We diagram the possible channels through which a part¬ 
nership can be terminated in Figure 5, and show the destruc¬ 
tive close encounters with the IMBH on the cluster evolution 
timeseries in Figure 6. 

In Section 4, we varied several key cluster parameters to 
consider the sensitivity of IMBH companions to cluster and 
IMBH parameters. We considered clusters with a higher stel¬ 
lar remnant retention fraction, smaller IMBH mass, and larger 
particle number as compared to our fiducial parameter set. 
These comparisons show varying compositions of compan¬ 
ion stellar type demographics (Figure 10). In particular, the 
retention of more NSs and BHs leads to a substantially larger 
fraction of BH companions as these objects mass-segregate 
to interaction with the IMBH (e.g. Leigh et al. 2014). These 
demographic differences propagate to minor changes in typi¬ 


cal companion orbital semi-major axis (or period) as shown in 
Figure 9, where we break up the semi-major axis distribution 
of companions by stellar type. We show that companion hier¬ 
archy (as defined by aoja\) depends on IMBH mass. Lower- 
mass IMBHs show companions that are more highly segre¬ 
gated from the remaining stellar system than do more massive 
IMBHs. This suggests that the capture of hierarchically seg¬ 
regated companions is a unique property of the IMBH mass 
range. 

The segregation of close companions to IMBHs, as a prop¬ 
erty unique to BHs in the intermediate mass range, has 
dramatic implications for the production of repeated fiaring 
episodes in tidal interactions between stars and IMBHs, as 
seen in the remarkable object HLX-1 (Farrell et al. 2009). 
A combination of stellar dynamical segregation, secular per¬ 
turbation from less bound stars, and tidal interaction with 
the IMBH (e.g. Hopman et al. 2004; Baumgardt et al. 2006; 
Blecha et al. 2006) could place a star in an orbit similar to the 
1 yr period observed for HLX-1. Continued interactions with 
the cluster background could then account for the perturbation 
of this orbit into a disruptive configuration. 

In closing, we note that we may be able to leverage detailed 
knowledge of IMBH companions to reveal or constrain the 
still-uncertain presence of IMBHs in dense clusters includ¬ 
ing GCs, young clusters, and compact galactic nuclear clus¬ 
ters. IMBHs themselves are notoriously tricky to uncover, but 
our analysis suggests that they may reveal themselves several 
transient channels. Gravitational waves generated by IMRI 
events should carry large enough amplitudes to be detectable 
to moderate redshifts. Our results indicate that these events 
occur, even in the modeled low-mass cluster, at rates similar to 
analytic predictions (Mandel et al. 2008), offering substantial 
promise for their detection by either ground or space based 
laser interferometers. Electromagnetic emission might arise 
through steady-state accretion from companion winds. Our 
calculations, shown in Figure 12, suggest that this may ex¬ 
ceed the IMBHs Bondi accretion rate from the cluster medium 
only a small fraction, ^ 5%, of the time, but it none-the-less 
suggests that a lower limit accretion rate exists for IMBHs in 
dense clusters. Finally, we discuss extremely luminous tran¬ 
sient accretion occurring in rare episodes of tidal disruption 
of stars by the IMBH. Especially if the IMBH launches jets 
in these events, these transients could trace the presence of 
IMBHs in clusters to cosmological distances. 
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